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Abstract 

A systematic study of the total /3-decay half-lives and /3-delayed neutron 
emission probabilities is performed. The /3-strength function is treated within 
the self-consistent density-functional + continuum- QRPA framework includ- 
ing the Gamow- Teller and first-forbidden transitions. The experimental total 
/3-decay half-lives for the Ni isotopes with A <76 are described satisfacto- 
rily. The half- lives predicted from A=70 up to ^4=86 reveal fairly regular 
yl-behaviour which results from simultaneous account for the Gamow- Teller 
and first-forbidden transitions. For Z w 28 nuclei, a suppression of the delayed 
neutron emission probability is found when the iV=50 neutron closed shell is 
crossed. The effect originates from the high-energy first-forbidden transitions 
to the states outside the Qp — ^-window in the daughter nuclei. 



1 



Typeset using REVTpjX 



I. INTRODUCTION 



Physics of the complex nuclear systems with high isospin asymmetry is a rich and multi- 
faceted field. It involves such exotic objects like neutron stars and short-lived /5-unstable 
nuclei with a large neutron (proton) excess. The ground state and /3-decay properties of these 
radioactive isotopes give a valuable information on structural evolution far off stability and 
also provide input for supernova explosion calculations. In particular for the modeling of the 
r-process nucleosynthesis, the most important are the nuclear masses which define the path 
of the r-process through the neutron-rich domain of the nuclear chart. The /3-decays and 
v{p) captures are also essential, as they regulate the flow of the material to high Z- values and 
set up the r-process time scale. The nuclei near the new doubly-magic 78 Ni, 132 Sn and "east" 
of 208 Pb are of special importance. Spectacular progress of the RNB experiments in these 
regions has been reached recently due to Z-selective resonance-ionization laser ion-source 
technique (see e.g [1]). It has provided a unique testing ground for theoretical approaches 
to exotic nuclei and modeling of explosive stellar events. 

In this short paper we perform the microscopic study of the total /5-decay half-lives and 
/5-delayed neutron emission probabilities of very neutron-rich nuclei. The simultaneous anal- 
ysis of these /3-decay observables gives a possibility to reconstruct the /3-strength function. 
Specifically, one may gain an insight about the relative contribution of its Gamow- Teller 
(GT) and first-forbidden (FF) components. A possibility of competition between the GT 
and FF decays in nuclei with the neutron-excess bigger than one major shell has been of- 
ten discussed in the literature (see [2] and Refs. therein). The experimental evidence of 
that has been known since the first measurements in the region of 132 Sn [3]. However, the 
microscopic global calculations [4,5] have treated the total /3-decay half-lives and delayed 
neutron emission probabilities (P n ) in allowed transitions approximation. (The estimated 
contribution of the unique FF transitions from [6] has to be taken with some reservations.) 

Clearly, the single /9-decay channel approximation is not adequate for description of the 
isotopic dependence of the /3-decay characteristics, especially for the nuclei crossing the 
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closed N and Z shells. It has been pointed out recently in [8] that self-consistent framework 
of the finite Fermi system theory including the GT and first-forbidden FF decays on the same 
footing allows for reasonably sound predictions of the ground state properties and /3-decay 
characteristics of very neutron-rich nuclei. An important contribution of the first-forbidden 
(FF) non-unique decays to the total half-lives of neutron-rich nuclei above the closed proton 
and neutron shells has been stressed in [8,9]. The purpose of the present work is to examine 
the impact of the high-energy FF decays on the /3-delayed neutron emission in the nuclei 
near the closed shells at Z=28 and N=50. In what follows we briefly outline the particulars 
of the theoretical framework, present new calculations of the total /3-decay half-lives and 
/5-delayed neutron emission probabilities and discuss possible uncertainties. 



II. THEORETICAL ANALYSIS 

The /^-delayed neutron emission is basically a multi-step process consisting of a) the 
/9-decay of the precursor (A, Z) which results in feeding the excited states of the emitter 
nucleus (A, Z + 1); b) neutron or/and 7-emission to an excited state or to the ground state 
of the final nucleus (A — 1, Z + 1). The difference in the characteristic time-scales of the f3- 
decay and subsequent particle emission processes justifies an assumption of their statistical 
independence. Thus, the total probability of the delayed neutron emission accompanying 
the /5-decay to the excited states in the daughter nucleus can be expressed as [10] 




where u is the energy of nuclear transitions between the ground state of the parent nucleus 
and final states of the emitter 1 . The neutron emission threshold of the emitter B n (Z + 1) 
and total /3-decay energy Qp are calculated self-consistently [7]. 

1 Within the self-consistent QRPA-approach one actually calculates the transition energy uo 
(clearly, an additional "outer" variable Qp [4] is needed to "derive" E x ). The transition ener- 
gies for the main GT and FF decays reveal a smooth A-dependence which simplifies the analysis. 
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In the Eq.l, Ti/ 2 stands for total (3-decay half-life 

1/T 1/2 = D- 1 (G A /G V ) 2 / du f (Z + 1,oj)J2 < kj > Spiu;^) . (2) 

Here the integrated lepton (e~, 9 e ) phase-space volume for allowed transitions (the inte- 
grated relativistic Fermi function) depending on the electron energy W 

fWmax 

fo= F(Z + 1, W) p W (W max — W) d W (3) 

J m e c 2 

takes into account for the Coulomb and finite-size corrections [11]. 

If the /3-decay transition energy relative to the mother nucleus ground state satisfies 
uj < Q/3 — B n = Qj3 n , the neutron decaying state locates within the so-called Qp n - window 
in the particle continuum; the energy of the emitted neutron being E n = Qp — B n — uj — E 7 . 
The particle emission probability from the level \j^) of the emitter to the level |jJ J ) of the 
final nucleus, as a function of the emitted neutron energy E n is given by 

p{r ^ En) = T n (If+l(E,y (4) 

where T n (j n i j, E n ) is the neutron escape width and T 7 (j 7T i j, E 7 ) - gamma ray width of the 
emitter. These can be expressed through the transmission coefficients for 1-wave neutron 
emission with the total angular momentum j n and corresponding 7-transmission coefficients 
calculated from statistical model of Hausser-Feshbach. In what follows we will assume 
that the 7-emission from neutron- unbound states is absent: T 7 <^ T n . The approximation 
simplifies the calculations but may cause some overestimation of the resulting P n -values. 

The intensity of feeding the excited states in the daughter nucleus is given by the set 
of /3-strength functions Sp{uj^) which describe the spectral distributions of the matrix 
elements of the GT and FF /3-decay transitions. The index (3 in the sum corresponds to the 
Gamow- Teller (GT) term with L = 0, J = 1 and non-unique first-forbidden (FF) terms with 
L — 1, J = 0, 1 treated in the ^-approximation. The unique J = 2 term can be also retained 
but in the ^-approximation it should be of minor importance. For the GT and non-unique 
FF decays < kj=o,i >=1, and for the unique FF decays, < kj =2 >—fi/fo, with the Fermi 
integral f\ calculated as in [12]. 



To calculate the /3-strength functions Sp(uj,^) we use the continuum QRPA approach 
based on the self-consistent ground state description within the local energy-density func- 
tional (DF) theory. The Fayans phenomenological DF [13] consisting of a normal and a 
pairing part is adopted. The DF3 version of the functional [7] contains the two-body spin- 
orbit and velocity dependent effective NN-interactions important for the full consistency, as 
well as the isovector spin-orbit force. The latter ensures a correct description of the single- 
particle levels near the "magic-cross" at 132 Sn [14]. We have studied also the ground state 
description given by the Skyrme MSk7 force [15]. The like-particle T = 1 pairing in the 
ground state is treated in the diagonal approximation on the basis confined by the E cuto ff =15 
MeV; the empirical strength of the T = 1 pairing decreases slightly with increasing A. 

The continuum QRPA (CQRPA) equations of the finite Fermi system theory are solved 
with an exact treatment of the particle-hole (ph) continuum. The width 7 includes both 
the escape width T 1 * for the decay to the continuum, as well as the spreading width of 
the isobaric states TK Within the CQRPA, the latter is assumed to depend linear on the 
excitation energy [16] (the cut-off at Au = 2.5T^ is employed for the high-energy wings of 
the individual excitations with u < Qp). 

A finite-range effective NN-interactions in the particle-hole (ph) channel is chosen in a 
5 + 7r + p form. The one-7r and one-p exchange terms modified by the nuclear medium 
are important in describing the magnetic properties of nuclei and the nuclear spin-isospin 
responses. The competition between the one-pion attraction g n Q < and contact spin- 
isospin repulsion g' Q > 0, determines a degree of "softness" of the pionic modes in nuclei 
that influences the /3-decay half-lives. The set of the NN-interaction parameters #'=0.98, 
(^=-1.38 has been used, corresponding to the quenching factor of the spin-isospin response 
function Q = e g [ar] 2 =0.81 (see [8]). 

The effective T = NN-interaction in the particle-particle (pp) channel is chosen in the 
same 5-function form, as the like-particle T = 1 pairing [8]. Neglecting the effective T = 
interaction would destroy the SO (8) symmetry of the QRPA equations and cause unrealistic 
odd-even staggering of total /3-decay half-lives. Importantly, the CQRPA-like equations of 
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the FFS allows for a reasonable description of nuclear spin-isospin modes with practically 
A-independent constant of the pp-interaction which is far from the instability point in the 
pp channel [17]. 

III. RESULTS 
A. The /3-decay half-lives 

The total /3-decay half-lives for Ni-isotopic chain are displayed at Fig.l in comparison 
to the finite-range droplet model [5] and experimental data [18] - [20]. The FRDM calcu- 
lations overestimate the experimental data and predict strong odd-even effect in the total 
/3-decay half-lives. This clearly demonstrates a significance of the effective pn-interaction 
in the particle-particle channel. Ignoring this interaction violates the SO(8)-symmetry of 
the QRPA equations and leads to unrealistic scale of the odd-even staggering of the total 
half-lives along the isotopic and isotonic chains. The improved calculation performed in the 
present work is in better agreement with the experimental data for Ni isotopes than the 
one from [8]. A slight underestimate of the experimental half-lives [18]- [20] is observed 
for A = 74 — 76. More importantly, the total half-lives in a rather long isotopic chain of 
A = 71 — 86 reveal a fairly regular behavior. Such a smooth A-dependence of the half-lives 
can be achieved only by taking into account the forbidden decays. To stress this point, we 
show at the Fig.2 the isotopic dependence of the energies for the main GT and FF transi- 
tions within the Q^-window. Evidently, the FF decays play a negligible role for the nuclei 
with A < 78 because of the low transition energy (small available phase-space). In contrast, 
after crossing the iV=50 shell, the high-energy FF transitions give a dominant contribution 
to the total half-life for the nuclei with A > 79. 

For the Ga isotopes with A < 84 (Fig. 3), both calculations give practically a similar 
agreement with the experimental data [21] - [23]. As seen from Fig.3, the FRDM+RPA 
calculations [5] underestimate the experimental data for A=84, while our calculations over- 
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estimate them. Apparently, this is related to the onset of the ground state deformation in 
the isotopes with A > 84 which has not been included in our calculations. 

B. The /3-delayed neutron emission 

Within the allowed transitions approximation used in the global calculations [4,5], the 
isotopic dependence of the P n -value is rather schematic. It is mainly defined by the the decay 
to the so-called GT pigmy- resonance located within the Q/3-window. The matrix element of 
the latter increase with N-Z, as well as the phase space Qp n available for the delayed neutron 
emission. Thus, the P n -value along the isotopic chain simply tends to the maximum after 
the GT pigmy- resonance turns to be located within the Q/3n-window. On the other hand, if 
the model includes the different /3-decay channels, an isotopic dependence of the P n -value is 
defined by the relative energies of the GT and FF transitions which behave differently with 
increasing A. Importantly, for the large enough neutron excess, the additional high energy 
FF transitions appear. 

The calculations of the /5-delayed neutron emission probabilities have been performed for 
several isotopic chains with Z ^28. Below, only the results for Ni and Ga isotopic chain are 
discussed in detail. As can be seen from Fig. 2, for the Ni isotopes with A <79, no neutrons 
can be emitted following the decay to the pygmy GT resonance and low transition energy 
FF states in the emitter nuclei. According to the present DF+CQRPA calculations they are 
located higher than the Q^n-window (Fig. 2). For this reason, the increase of total delayed 
neutron emission probability in the isotopes with A<79 (Fig.4) is entirely due to relatively 
low-energy GT and FF /5-decays. 

Starting from A=79, the GT pigmy-resonance enters the Q/3n-window (Fig. 2). Assuming 
a pure GT-decay, the P n -values for A >79 would tend to 100%. However, for the isotopes 
with A >79 crossing the iV=50 shell, the drastic change in the relevant shell configurations 
takes place. From Fig. 2 which shows the energies of the different FF transition with the 
total transferred momentum J = 0, it is seen that with filling of the z/2rf 5 / 2 -orbital in the 
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A >79 isotopes, the high-energy FF component of the /^-strength function appears due to 
the v2d b / 2 — > vr2/ 5 / 2 transition. As mentioned above, the high-energy FF transitions give 
a dominant contribution to the total /3-decay half-lives in the Ni isotopes with A >79. At 
the same time, their calculated energies are higher than the Q /^-threshold. Hence, the 
high-energy part of the /9-strength function does not affect the delayed neutron emission, 
the feature that translates into a "gap"-like pattern in the P„(A)-curve (Fig. 4). For the 
nuclei with A >86 the main J = transitions are located in the Q/3 n -window (Fig.2), and 
the corresponding P n - values tend to the maximum (Fig.4). It would be important to check 
our prediction for Ni isotopic chain, though for A >80 for which the difference with the GT 
approximation is observed it may be difficult at the moment. 

In the context of the possible experimental measurements, the isotopes with higher Z are 
of special interest. To exemplify this we consider the case of Ga isotopes. For 80-81 Ga, the FF 
transitions are of minor importance but with increasing A their contribution becomes very 
significant. It follows from Fig. 5 that in contrast with Ni isotopes, the P n -values calculated 
for 82-87 Ga isotopes are suppressed by the factor of 4-5 compared to the ones corresponding 
to the GT approximation. Notice that in 82 ~ 84 Ga, the evaluated P n -values [23] obtained 
by averaging of the existing experimental data [22] - [24] can be described within the GT 
approximation alone. However, it has to be realized that the discrepancy of the experimental 
data is very high. Interestingly enough, in 83 Ga, the ones from [24] are in agreement with the 
GT+FF calculations. Assuming that reliability of the experimental data for Ga isotopes has 
been questioned in [23], it would be of great importance to perform the new measurements 
for this isotopic chain. 

To conclude, one has to mention that in Z=28 region, the delayed neutron emitter nuclei 
have relatively high Qp- values. For that reason, the accuracy of the calculated Ti/ 2 and Pn- 
values is generally higher than for the heavier nuclei with lower Qp- values. However, as the 
delayed neutron emission is typically a threshold phenomenon, it is of interest to analyze the 
different factors which may influence its probability. The first concerns the nuclear pairing, 
as it changes the one-quasiparticle level energies and calculated P n -values. This issue has 
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been tackled in [26], however the only calculation was performed for 137 I with a very low 
Ql3 n — 1.6 MeV in which case a sensitivity to the paring strength is extremely high. It would 
be of importance to study in detail how the predicted effect might be influenced by the 
different prescriptions of nuclear pairing. To roughly estimate the possible sensitivity we 
varied the _B n -values within the ±1 MeV. It turns out that for the nuclei with high Qp n - 
values, the resulting impact factor is not that strong e.g. for 83 Ni the P n -value varies within 
the margins of 70.8 to 65.3% instead of 62.7% shown at Fig.4. 

Second, the "gap-like" behavior of the P n - values turns to be fairly robust against the 
allowed choice of the effective NN-interaction parameters. (Using the set of the parameters 
corresponding to the quenching factor of e g =0.8 instead of 0.9 [8] leads to the P n -value of 
71% for the 80 Ni.) 

Third, it is worth to estimate a possible influence of the /^-strength redistribution due 
to the deformation and/or the specific np-nh correlations which may shift the /5-strength 
outside or inside the Q/jn-window. As it is seen from Fig. 3 for the nuclei of interest beyond 
80 Ni, the GT pigmy-resonance energy is much lower than the position of the Q^-window and 
such a shift seems quite unlikely. As for possible impact of the deformation, the calculations 
[15] predicts the quasi-spherical shape ((3 2 < 0.1) for Ni isotopes with A = 74 — 94 and for 
Ga isotopes with A = 79 — 83. 

Finally, as a proper account for the particle-emission channels within the Hausser- 
Feshbach framework has not been performed, the predicted P„-values may be considered 
as an upper limit estimate. Moreover, for very neutron-rich nuclei the P n -values become 
low, as well as the density of the neutron unbound levels in the near threshold region. In such 
a case, a more detailed study of the /^-delayed neutron emission may be of need considering 
the contribution of direct neutron decay of the isobaric states. 
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IV. SUMMARY 



The DF+CQRPA model of the delayed neutron emission is elaborated with the Gamow- 
Teller and first-forbidden /3-decay modes taken into account. The calculations of the total 
/9-decay half-lives and delayed neutron emission probabilities has been performed for nuclei 
near the closed shells at Z=28 and N=50. Within the extended model, an agreement with 
the experimental data on the total half-lives is better than in our previous calculations 
[8]. A suppression is predicted of the /^-delayed neutron emission probabilities for Z ^28 
isotopes above the A=50 shell. The first-forbidden decays in these nuclei are found to give a 
dominant contribution to the total half-lives. A threshold character of the delayed neutron 
emission makes it sensitive to deblocking of the shell-configurations responsible for the high- 
energy forbidden /5-decay transitions. Inclusion of the FF /3-decays leads to reduction of 
the total delayed neutron emission probabilities, as the calculated energies of the main FF 
transitions are systematically higher than the Qp n - values. 

For Ni isotopes, the scale of predicted reduction of the P n -values (compared to the ones 
estimated in the GT approximation) is at the edge of the experimental precision achievable 
currently with the isotope-separated beams. A similar "gap-like" behavior of the total prob- 
ability of delayed neutron emission has been found in our calculations for the Cu isotopes. 
A stronger suppression of the P„-values predicted for Ga isotopes, as well as for the Zn and 
Ge-As isotopic chains provides a certain challenge, as these isotopes are more accessible for 
the experimental measurements in the short term [25]. 
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FIGURES 

FIG. 1. The total /3-decay half-lives for the Ni-isotopes predicted from the FRDM+RPA [5] 
and DF+CQRPA [8]; the experimental data are taken from [18] - [20] 

FIG. 2. The calculated position of the Q/3 n -wmdow for the delayed neutron emission compared 
to the energies of the GT pigmy-resonance and main J 11 = 0~ transitions. For the lowest energy 
branch of the 0~ transitions (a>i), the points for A=77-81 are shown only 

FIG. 3. The total /3-decay half-lives for the Ga-isotopes predicted from the FRDM+RPA [5] 
and DF+CQRPA ( present work ); the experimental data are taken from [23] 

FIG. 4. The total /3-delayed neutron emission probabilities for the Ni isotopic chain calculated 
within the DF+CQRPA for the allowed GT transitions only (A) and in the GT+FF approximation 

(O) 

FIG. 5. The total /3-delayed neutron emission probabilities for the Ga isotopic chain calculated 
within the DF+CQRPA for the allowed GT transitions only (A) and in the GT+FF approximation 
(O)- The experimental data are: a) [23], b) [22], c) [24] 
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